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Abstract 

In this paper we derive the exact solution of the multi-period portfolio choice problem 
for an exponential utility function under return predictability. It is assumed that the asset 
returns depend on predictable variables and that the joint random process of the asset 
returns and the predictable variables follow a vector autoregressive process. We prove 
that the optimal portfolio weights depend on the covariance matrices of the next two 
periods and the conditional mean vector of the next period. The case without predictable 
variables and the case of independent asset returns are partial cases of our solution. 
Furthermore, we provide an empirical study where the cumulative empirical distribution 
function of the investor's wealth is calculated using the exact solution. It is compared 
with the investment strategy obtained under the additional assumption that the asset 
returns are independently distributed. 
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1 Introduction 



Investment analysis and portfolio choice theory are very important and challenging topics in 
finance and economics nowadays. Since Harry Markowitz (1952) presented his mean-variance 
paradigm portfolio theory has become a fundamental tool for understanding the interactions 
of systematic risk and reward. 

It is well known that the mean- variance optimization problem of Markowitz (1952) is equiv- 
alent to the expected exponential utility optimization under the normality assumption (see 
Merton (1969)). Unfortunately, his approach only gives an answer to the one-period portfolio 
selection problem in discrete time but it says nothing about the multi-period (long-run) setting. 
Therefore, it is of importance to investigate the multi-period portfolio optimization problem 
which is of great relevance for an investor as well. The multi-period portfolio selection problem 
has been analyzed for different types of utility functions by, e.g., Mossin (1968), Merton (1969, 
1972), Samuelson (1969), Elton and Gruber (1974), Brandt and Santa-Clara (2006), Basak and 
Chabakauri (2010). 

The continuous case has already been solved for many types of utility functions in the 
one- and multi-period case by Merton (1969). A number of generalizations under weaker 
assumptions are given, among others, by Ait-Sahalia et al. (2009) and Skaf and Boyd (2009). 
Exact solutions in discrete time are even more difficult to obtain in the multi-period case. 
Mossin (1968) considers the case of one risk and one riskless asset. He derives conditions 
when the multi-period strategy becomes myopic or partial myopic. Frequently this can be 
achieved by demanding independent asset returns. However, the assumption of independence 
is unfortunately not fulfilled in many applications. 

For an investor it would be very useful to have a closed-form solution of the discrete multi- 
period portfolio problem. Moreover, it is desirable that the optimal portfolio weights can 
easily calculated at each period. Analytical solutions of the multi-period optimal portfolio 
choice problems are hard to obtain and they are only available for some partial cases which 
are often derived under very restrictive assumptions on the distribution of the asset returns. 
For instance, a closed-form solution for the multi-period portfolio choice problem exists for 
the quadratic utility function under the assumption that the asset returns are independently 
distributed (see Li and Ng (2000), Leippold et al. (2004)). 

In the present paper we consider an investor who invests into k risky assets and one riskless 
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asset with an investment strategy based on the exponential utility function 

U(W t ) = -e~ aWt . (1) 

Here W t denotes the investor's wealth at period t and a > stands for the coefficient of 
absolute risk aversion (ARA), which is a constant over time for the exponential utility (CARA 
utility). The application of the exponential utility function is more plausible than the use of the 
quadratic utility since the first one is monotonically decreasing. That is why the exponential 
utility function is commonly used in portfolio selection theory. Moreover, the optimization 
of the expected exponential utility function leads to the well known mean-variance utility 
maximization problem and consequently its solution lays on the mean- variance efficient frontier. 

We derive a closed-form solution of the multi-period portfolio choice problem with the ex- 
ponential utility function ([I]) under the assumption that the asset returns depend on certain 
predictable variables. The joint process consists of the asset returns and the predictable vari- 
ables and it is assumed to follow a vector autoregressive (VAR) process. This approach is very 
popular in finance and it is often used for modeling the asset returns (see, e.g., Campbell (1991, 
1996), Barberis (2000), Brandt (2010)). 

The rest of the paper is organized as follows. In Section 2, the main result of the paper 
is given. In Theorem 1 an analytical expression of the portfolio weights is provided for each 
period. In Corollary 1, the case without a predictable vector is considered while independent 
asset returns are treated in Corollary 2. In Section 3 a short empirical study is presented. The 
performance of the derived strategy is compared with the one for independent asset returns. 
The comparison is performed in terms of the cumulative empirical distribution function of the 
investor's wealth at the end of the investment period. We find significant improvements if the 
dependence structure is taken into account. Section 4 contains a short summary. 



2 Mult i- Period Portfolio Problem for an Exponential 
Utility 

In this section we derive the analytical solution of the multi-period portfolio choice problem for 
an exponential utility function assuming that the asset returns and the predictable variables 
follow a VAR(l) process. 

There are only a few papers in literature where the exponential utility function is considered 
in the multi-period discrete time setting. For instance, Canakoglu and Ozekici (2009) solved the 
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portfolio choice problem assuming that the stochastic market follows a discrete time Markov 
chain and all parameters of the asset returns, i.e., mean vector and covariance matrix, depend 
only on the current state of the stochastic market and not on the previous states which is 
equivalent to the assumption of independence in our settings. In the paper of Soyer and 
Tanyeri (2006) a Bayesian computational approach with the exponential utility was presented. 
The authors write that the solution of the multi-period portfolio choice problem with the 
exponential utility under the assumption of normality "... cannot be evaluated in closed form 
and the optimal portfolio cannot be obtained analytically" . In this paper, however, we present 
an exact solution assuming that the asset returns follow a vector autoregressive process with 
predictable variables. 

Let X( = (X t> i, X tj 2, ■ ■ ■ , Xtje) denote the vector of the returns of k risky assets and let r/^ 
be the return of the riskless asset at time t. Let z t be a p-dimensional vector of predictable 
variables. We assume that Y t = (X.' t ,z' t )' follows a VAR(l) process given by 



with e t ~ A/"(0, where is a positive definite deterministic matrix function. Let 

J~t denote the information set available at time t. Then Y t \J-t-i ~ J\fk +p (p, t ,'S(t)), i.e., the 
conditional distribution of Yt given J-t-i is a k + p dimensional normal distribution with mean 
vector fi t = E^^Ft-i) = E t _i(Y t ) and covariance matrix Var(Y t | J-'t-i) = 

The stochastic model ^ is described in detail by Campbell et al. (2003) who argued that 
the application of VAR(l) is not a restrictive assumption because every vector autoregression 
can be presented as a VAR(l) process through an expansion of the vector of state (predictable) 
variables. The idea behind this approach is to find a vector of predictable variables z t which 
is mostly correlated with the asset returns and to build a VAR(l) process with respect to the 
asset returns X t and the vector of predictable variables z t . The choice of z t depends rather on 
the data and not on the utility function. Possible predictable variables are, e.g., the dividend 
yield (cf. Campbell at al. (2003)), the term spread (see, e.g, Brandt et al. (2006)) or another 
asset return. 

From ([2l we obtain the following model for X t expressed as 



where is a k x k identity matrix and Ok, P is a k x p matrix of zeros. Consequently, X t | J-'t-i ~ 



A4(Mt, £(*)), where fi t = E{X t \J r t -i) = v + ^Y t ^ and £(t) = Var(X t |JF t _ 1 ) = LE(t)I/. 



Y t = v + $Y 4 _! + e t 



(2) 



Xt = Lj> + L$Y t ! + Lit = v + $Y t _! + e t with L = [I fc O k)P ) , 



(3) 
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Let w t = (wt t i, w tj 2, ■ ■ ■ , Wt,k)' denote the vector of the portfolio weights of the k risky assets 
at period t. Then the evolution of the investor's wealth is expressed as 

W t = W t -! (1 + r u + wU(X t - r u l)) = W t -i (Rf,t + <-!%) , (4) 

where Rf t t = 1 + r/^ and X t = X t — with fl t = £? t _j(Xt) = f + 3>Y t _i — The aim 

of the investor is to maximize the expected utility of the final wealth. 
The optimization problem is given by 



^(0,^0,^0)= max E t [U{W T )} (5) 
with the terminal condition 

U(W T ) = -exp(-aW T ) for a > 0. (6) 

Following Pennacchi (2008) the optimization problem ^ can be solved by applying the follow- 
ing Bellman equation at time point T — t 

V(T -t.WT-t^T-t) = max£ T _ t [ max E T _ t+1 [U(W T )}\ (7) 

WT-t 



{w s } T " 



s=T-t+l 



max Ex-t 

WT-t 



V(T -t+1, W T -t (rf,T-t + w T-t+i x T-t+i) , Fr-t+i 

subject to ([6]), where w^_ t+1 are the optimal portfolio weights at period T — t + 1. Note that 
in contrast to the static case now the vector of optimal portfolio weights Wx-t+! is a function 
of the weights of the next periods, i.e., of w T _ t+1 ,w T _ t+2 , . . . ,Wt-i, what is the consequence 
of the backward recursion method (see, e.g. Pennacchi (2008)). 
For the period T — 1 we get 

V(T-1,Wt-!,T T -!) 



exp(— aWT-i(Rf,r + w^jXy)) 



= — exp(—aWT-iRf,T)ET-! [exp(— qWt-iw^jXt)] 

= exp(— aWx~!Rf,T) exp —a(WT—iyfr-it^T — — •w' t _ 1 'S(T)-wt-iW^_ 1 ) ^ — > max. (8) 
The last optimization problem is equivalent to 

WT-iWr-ifi'T ~ -^ w 't-i^(T) w t-iWt-! — > max over w^i . (9) 

Taking the derivative and solving ^ with respect to w^-i we get the classical solution for the 
period T — 1 

w* T _, = -^—^-\T)fl T = —^—(Lt^LT'^T + ^YT^) with i> T = v-r f)T l. (10) 
In Theorem 1 the multi-period portfolio weights for all periods from to T — 1 are given. 



Theorem 1. Let X r = (X T) i,X Tt 2,...,X T} k)' be a random return vector of k risky assets. 
Suppose that X T and the vector of p predictable variables z T jointly follow a VAR(l) process 
as defined in M). Let r/ iT be the return of the riskless asset. Then the optimal multi-period 



portfolio weights are given by (10) for period T — 1, 



Wt-2 



->)RfT V 



aW / r_2-R/,r 



(11) 



and 



w 



T-t — T 

aW T _ t n Rf,i 

i=T-t+2 



LS X {T -t + l)fi* T _ t+1 -L*'S 1 (T-< + 2)(^_ t+3 +r /iT -t+2*L'i; 



(12) 



OTi/i Ar-t+i = Ar-t+i - r /iT _ t+2 L'l and i>T-t+3 = ^ - r/.T-t+sL'l , 



fort = 3,...,T. 



Proof. The value function at time point T — 2 is obtained by substituting (|10|) into (|8 
V(T-2,^ T _ 2 ,J- T _ 2 ) 



T-2 



exp ( -aW T -\Rf,T ~ -St 



exp (— aW / r-2-R/,T-i-R/,r) -Er- 



exp ( — aWr-2-R/,TWy_ 2 XT-i — —St 



(13) 



where = jfi' T Ti{T) 1 (i T . Second, according to the properties of VAR(l) processes we get 
that 

s T = Y' T _ 1 &?:(T)- 1 ®YT-i + 2Y' T _ 1 &?:(T)- 1 i> T + , 

(14) 

with Y t = Y t — r/^L'l. This is a quadratic form with respect to the conditional normally 
distributed vector Y T _ 1 . Moreover, using LL' = I fc , 



Wt oX 



T-2 A T-1 — w T-2l- Ljl T-l 



fLY 7 



r/,rl) 



w 



(15) 



and (13), we get 



V(T - 2, W^-s, J" T - 2 ) = - exp (-aW T _ 2 Rf,TRf,T- 



x 



T-2 



exp ( —^Y' T _ 1 B(T)Y T _ 1 - b(w T _ 2 )'Y T _ 1 - c 



(16) 



where B T = ^(T)" 1 *, b(w T _ 2 ) = ^'^(T^Pt + aW T - 2 Rf,T^T-2 +r f>T B(T)L'l and 
c = i(i> T + r/ ) r*L'l)'E(T)- 1 (i/ r + r /iT $LT). 



Following Mathai and Provost (1992, Theorem 3.2a.l) the expectation given in (16) is the 
moment generating function of the quadratic form in normal variables at point —1. Hence, it 
holds that 

V(T-2,W T _ 2 , JF r _ 2 ) = -expi-aWT^Rf^Rf^) |I + B(T)E(T - 1)|"3 
1 



x exp 



2 (At-i " r /iT L'l)'S(T - IJ-^At-x " r f , T Vl) - c 



+ g (Ar-i - r f)T L'l - E(T - l)b(w T - 2 ))' 

x (I + B(T)S(T - l))" 1 !] - '(T - 1) (At-i " r/.rL'l - E(T - l)b(w r _ 2 ))] , (17) 

where E t _ 2 [^t-i] = At-i — r /,TL'l. Thus, the optimization problem V(T — 2, Wt-2,3~t-z) — > 
max is equivalent to V*(T — 2, Wt-2, Ft -2) ~ ► max, where 

V*(T - 2, W T _ 2 , J- r _ 2 ) = -^(At-x - r f , T L'l - E(T - l)b(w T _ 2 ))' 

(I + B(T)E(T - l))- 1 ^^ - l)- 1 (/i^ - r f>T L'l - E(T - l)b(w r _ 2 )) 

= -\(fcr-i ~ r f , T Vl - E(T - l)b(w T - 2 ))' 

x (t(T - 1) + E(T - 1)B(T)E(T - 1)) ' (At-i - r/.rL'l - E(T - l)b(w T - 2 )) • (18) 

Because the matrix E(T — 1) + S(T — 1)B(T)E(T — 1) is positive definite the maximum of 
V*(T — 2, Wt-2, J~t-2) is attained at w^_ 2 for which 

fi T _ x - r f)T L'l = £ T -ib{w T _ 2 ). (19) 

Using that LL' = I*. we obtain 

-i— (LS~ 1 (T-l)(AT-i-r/, T L'l)-L* / E- 1 (T)(i> r + r /)r *L'l)) . (20) 

T—2^f,T ^ ' 



T-2 



Furthermore, the equality (19) leads to 

U(T-2,H/ T _ 2 , Jr_ 2 ) = -|I + B(T)S(T-l)|-3exp(-c) 

x exp ( -aW T - 2 %r-ii2/,T - ^(Ar-i " r /ir L'l)'Xf X (r - l)(Ar_i - r/, T L'l)) (21) 
= -|I + B(T)E(T - 1)|~3 exp (-c) 

x exp f-H' T _ 2 B(T - 1)Y T _ 2 - b'(w T _ 3 )Y T _ 2 - c) , (22) 

where B(T-l) = ^'^(T-l)- 1 *, b(w r _ 3 ) = &t(T-iy 1 (i>-r ftT L , l)+aW T - 3 R fiT - 1 L , w T ^ 3 + 
r /)T -iB(T- l)L'l and c = \{u - r /)T L'l +r /iT _ 1 $L'l) / E(T - l) _1 (i> ~ r/,rL'l + r- /)T _i$L'l). 



Taking the conditional expectation from the value function (21 ) with respect to J-'ts we receive 



V(T - 3, W T - 3 , Tt-z) = -|I + B(T)S(T - 1)|~3 exp (-c) 

x £ T _ 3 (exp (~Y^_ 2 B(T - 1)Y T _ 2 - b'(w T -3)Y T -i - cj J . (23) 

Consequently, the value function for the period T — 3 has a similar structure than V^(T — 
2, Wt-2, J~t-2) with the only difference that B(T) is replaced by B(T— 1), b(w T _ 2 ) by b(w r _ 3 ) 
and c by c. It follows immediately that the optimal portfolio weights at period T — 3 are 

W ^" 3 = aW T ^R ST R fT ^ ( Lir V - 2)(M T _ 2 - r /)T -iI/l) - L*'E _1 (T - + ry.T-jSL'l)) (24) 

with i/y = v — r^yl/l. 

The last step is to use mathematical induction with basis T — 3 in order to receive the 
statement of Theorem 1. □ 

The results of Theorem 1 show us that the optimal portfolio weights at every period of 
time except the last one depend on the covariance matrices of the next two periods and the 
conditional mean vector of the next period. This property turns out to be very useful if we 
want to calculate the optimal portfolio weights for a real data set. 

Note that the case without predictable variables is a special case of Theorem 1. In this case 
the following expressions are obtained. 

Corollary 1. Let X r = (X T)1 , X r 2 , . . . , X r ^)' be a random return vector of k risky assets which 
follows a VAR(l) process as defined in but without a vector of predictable variables z T . Let 
r/ )T be the return of the riskless asset. Then the optimal multi-period portfolio weights for period 
T — 1 are given by 

1 1 

S" 1 (T)/i r = — — tT (T)(i/ T + *Y T _i) with i> T = u - r fjT l (25) 



w 



T 1 aW T -i K /r ~ J olWt-i 
and for t = 2, . . . , T by 

w*T- t = \ {^\T - t + l)fi T _ t+l - *'TT\T - t + 2)(i> T _ t+2 + r /jT - f + 2 $l)) 



aW T _ t fl R u 

i=T-t+2 



(26) 



Proof. The results of Corollary 1 are obtained in the same way as the results of Theorem 1 by 
putting L = 1^ and Y t = X 4 {p = 0) . □ 

In Corollary 2 the return vectors are assumed to be independent. 



Corollary 2. Let X T = (X T) i, X Tj 2, . . . , X T ^)' be a sequence of the independently and identically 
normally distributed vectors of k risky assets, i.e., X r ~ A/"(/x, S). Let r/ )T be the return of the 
riskless asset. We assume that S is positive definite. Then for all t = 1, . . . ,T the optimal 
multi-period portfolio weights for period T — t are given by 

Wt- ( = E -1 /* with jx = /j,- r />T _ i+2 l • (27) 

aW T - t n Rf,i 

i=T-t+2 

Proof. Corollary 2 immediately follows from Corollary 1 putting <3? = 0, E(t) = E and v = 
H- □ 

The results of Corollary 2 can be obtained as a partial case of Canakoglu and Ozekici (2009), 
where the stochastic market was presented by a discrete time Markov chain. In that case the 
asset returns depend on the present state of the market and not on the previous ones which 
implies the independence of the asset return over time. 

It is noted that the dynamics of the optimal portfolio weights in Corollary 2 is hidden in 
the coefficient of the absolute risk aversion a which is given by a T = I aWr- T ]~I Rf,i ) 

V i=T-r+2 ' J 

Moreover, the expressions of the weights themselves are proportional to the weights of the 
so-called tangency portfolio (cf. Ingersoll (1987, p. 89), Britten- Jones (1999)). Because the 
tangency portfolio is usually considered as a market portfolio in the single-period allocation 



problem (see, e.g., Britten- Jones (1999)) we treat the weights given in (27) as the weights of a 



benchmark portfolio in our empirical study presented in the next section. 



3 Empirical Study 

In this section we apply the results of Section 2 to real data. In following we consider an investor 
who invests into an international portfolio. The portfolio consists of the capital market indices 
of five developed stock markets, namely Belgium, Germany, Japan, the UK, and the USA. 
We deal with weekly data of the MSCI (Morgan Stanley Capital International) indices for the 
equity market returns from the January 4,2002 to December 4,2009. A process is fitted to the 
return series 

Xj = v + $Y t _x + e t with e t ~ UAT{0, S e ) . (28) 
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We get 
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-0.0889 


0.00988 


-0.441 


0.382 


2.79e 


-05 




0.0766 


-0.0643 


-0.03049 


-0.114 


0.133 



, and (29) 



0.0013085186 0.0010544496 0.0004365753 0.0009120373 0.0006781289 
0.0010544496 0.0013833540 0.0005648237 0.0010218539 0.0008332314 
0.0004365753 0.0005648237 0.0007994341 0.0004733366 0.0003667012 
0.0009120373 0.0010218539 0.0004733366 0.0010176793 0.0006927251 
0.0006781289 0.0008332314 0.0003667012 0.0006927251 0.0007242233 
It is remarkable that the last column of the matrix <fr has the largest values which indicate 
on a strong positive correlation between the US market and the other markets. Moreover, it 
shows that the influence of the US market on the return indices is larger than those of the 
domestic ones. Following Campbell et al. (2003) we choose the stock index of the US market 
predictable variable z t in our empirical study. 
Next, we calculate the weights of the two multi-period portfolio strategies for an exponential 
utility function. We want to compare the case of correlated return vectors given in Theorem 1 
with the case of independent variables given in Corollary 2 which completely ignores the time 
dependence structure well documented for real data. 

The performance of both strategies is compared with each other via an extensive simu- 
lation study based on 10 5 independent repetitions. The multi-period portfolio strategies are 
constructed for T e {13,26,52,104} and for the coefficient of relative risk aversion (RRA) 
a r G {0.8,2}. The RRA a r = aW is chosen as a constant absolute risk aversion (ARA) a in 
this study (without loss of generality we put Wo = 1)- In order to compare the performance 
of these two strategies we determine the empirical cumulative distribution function (ECDF) of 
the investor's terminal wealth for each strategy. 

The obtained results are presented in Figure 1 and 2. If we compare the performance of two 
portfolio strategies by their ECDFs, we should choose the strategy whose distribution function 
lies below the other because the probability of getting a larger wealth is larger for the strategy 
with a stochastically smaller distribution function. The strategy based on the weights given 
in Theorem 1 is denoted by EXP, while the notation EXP-iid is used for the method with the 
weights of Corollary 2. 
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Figure 1 presents the results for a smaller value of the coefficient of the relative risk aversion. 
We observe that EXP overperforms EXP-iid for all considered investment periods T. For 
instance, for T = 104 the probability of getting a wealth between 60 and 80 is equal to roughly 
25% for EXP while it is almost zero for EXP-iid. For a small horizon T there exists a small 
probability of bankruptcy for both strategies but it differs not significantly. For T > 52 the 
probability of a loss tends to zero. The EXP and the EXP-iid strategies both improve as T 
becomes larger what indicates their good performance in the long-run setting. 

Similar results are obtained for larger values of the coefficient of the relative risk aversion 
(see Figure 2). The performance of the EXP strategy is better for all T. From the other 
side, the probability of obtaining a larger value of the wealth is for both strategies smaller in 
comparison to the results presented in Figure 1. 

Using the results of both figures we can conclude that the EXP strategy has a higher perfor- 
mance for all T and risk levels. Of course this is not surprising since more information about 
the distribution of the asset returns is taken into account. On the other hand, ignoring the 
time dependence of the asset returns weakens the results with respect to the final wealth but 
it does not influence the probability of being bankrupt at the end of the investment period. 
Moreover, it has to be noted that the comparison of the ECDFs of the expected utilities is not 
relevant in our study because both strategies give the maximum expected utility in most of the 
cases and do not differ significantly. 

4 Summary 

Although the first formulation of the multi-period portfolio choice problem was already provided 
by Markowitz (1952), there are only a few results on closed-form solutions available in litera- 
ture. They are mostly derived under the assumption that the asset returns are independently 
distributed. Merton (1969) discovered that the maximization of the exponential utility function 
for normally distributed returns is equivalent to the maximization of the mean-variance utility 
function. Canakoglu and Ozekici (2009) obtained a closed-form solution for the exponential 
utility function under the assumption that the asset returns are independent. In general, the 
derivation of an analytical solution of the multi-period portfolio choice problem with the expo- 
nential utility for discrete time was considered to be very difficult (see, e.g., Soyer and Tanyeri 
(2006)). 

In the present paper we derive an exact solution of the multi-period portfolio selection 
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problem for an exponential utility function which is obtained under the assumption that the 
asset returns and the vector of predictable variables follow a vector autoregressive process 
of order I. Under the assumption of independence the obtained expressions of the weights 
are proportional to the weights of the tangency portfolio obtained as a solution in the case 
of a single-period optimization problem. We show that only the coefficient of absolute risk 
aversion depends on the dynamics of the asset returns in this case. The weights of the optimal 
portfolio derived without a vector of predictable variables are obtained as a partial case of the 
suggested general solution. In an empirical study we compare the derived multi-period portfolio 
strategies for real data taking five developed stock market indices. A very good performance 
of the general solution is observed which always overperforms the weights derived under the 
additional assumption that the asset returns are independent. 

The obtained results can be further extended by taking into account the uncertainties about 
the parameters of the data generating process. The analytical expressions of the weights can 
be used to derive the expected mean vector and the covariance matrix of the estimated weights 
which provide us the starting point for the detailed analysis of their distributional properties. 
This problem is not treated in the present paper and it is left for future research. 
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